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Dominant decay of a SM-like Higgs boson into particles beyond those contained in the minimal 
supersymmetric standard model has been identified as a natural scenario to avoid fine tuning in 
electroweak symmetry breaking while satisfying all LEP limits. In the simplest such extension, 
the next-to-minimal supersymmetric model, the lightest CP-even Higgs boson can decay into two 
pseudoscalars. In the scenario with least fine tuning the lightest CP-even Higgs boson has mass of 



\^ . order 100 GeV. In order to escape LEP limits it must decay to a pair of the lightest CP-odd Higgs 

o , 

\ bosons with Br(/i aa) > .7 and < 2mf, (so that a t^t or light quarks and gluons). The 

mass of the lightest CP-odd Higgs boson is controlled by the soft-trilinear couplings, Ax{mz) and 
^ • Af^{mz)- We identify the region of parameter space where this situation occurs and discuss how 

natural this scenario is. It turns out that in order to achieve ttIq < 2mb with Ax{mz), A^{mz) of 
^ ' order the typical radiative corrections, the required tuning of trilinear couplings need not be larger 

than 5-10 %. Further, the necessity for this tuning can be eliminated in specific SUSY breaking 
scenarios. Quite interestingly, Br(/i — s- aa) is typically above 70 % in this region of parameter space 
and thus an appropriately large value requires no additional tuning. 
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I. INTRODUCTION 



The Minimal Supersymmetric Standard Model (MSSM) is a widely studied possibility 
for physics beyond the standard model (SM). Its content in the matter and gauge sectors 
is fixed by all known particles and assumed superpartners. However, the choice of a two- 
doublet Higgs sector is made purely on the basis of minimality arguments. It is exactly the 
Higgs sector, namely the non-observation of the Higgs boson, that casts a shadow on the 
whole MSSM. It is becoming obvious that if the MSSM is the correct description of nature 
then the supersymmetric (SUSY) spectrum has to be quite unusual: heavy enough that 
sparticles escape direct detection, but not so heavy that electroweak symmetry breaking 
becomes unnatural. And, if the sparticles are not particularly heavy, the MSSM parameters 
must be special enough that the Higgs boson mass is pushed (through radiative corrections) 
above the experimental bound. ^ Another possibility, which does not require any special 
assumptions about the SUSY spectrum and parameters, is to abandon the minimal Higgs 
sector, for which there was never any deep reason anyway and which gives rise to the famous 
/i-problem. In this case, the expectations for Higgs phenomenology are modified and the 
tension from not having observed the Higgs boson at LEP can be eliminated ^. 

A particularly appealing extension of the SM or MSSM is the introduction of a completely 
new sector of particles which are singlets under the SM gauge symmetry. As such, this extra 
(E) sector would not spoil any of the virtues of the MSSM, including the possibility of gauge 
coupling unification and matter particles fitting into complete GUT multiplets. In addition, 
E-sector particles would have easily escaped direct detection. Of course, if this E-sector is 
completely decoupled from the SM then it plays no role in particle physics phenomenology 
at accelerators. However, it is possible that this sector couples to the MSSM through the 
Higgs fields. For example, the superpotential can contain a term in which the two Higgs 
doublets are combined in a SM-singlet form. Without additional Higgs fields, the coefficient 
of this form must have dimensions of mass, and the resulting superpotential component is 
the so-called yU.-term. In contrast, the E-sector can couple to this SM-singlet form in many 



^ Scenarios that lead to such a special SUSY spectrum were recently found, see for example mixed anomaly- 
modulus mediation [l| or gauge mediation with gauge messengers Q . Both scenarios generate large mixing 
in the stop sector which maximizes the Higgs mass, allowing all experimental limits to be satisfied with a 
fairly light SUSY spectrum. 
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ways, including at the renormalizable (dimensionless coupling) level. Such couplings would 
have a negligible effect on the phenomenology involving SM matter fields, whereas they can 
dramatically alter Higgs physics. For example, they would allow the lightest CP-even Higgs 
boson h to decay into two of these E-fields if the E-fields are light enough. In that case, the 
usual Higgs decay mode, h — > 66, might no longer be dominant and, since the hhh coupling 
is small, Br(/i bb) can be suppressed by a large factor. The strategy for Higgs discovery 
would then depend on the way the E-fields appearing in the decays of the h themselves decay. 
They might decay predominantly into other stable E-fields, in which case the MSSM-like 
h decays mainly invisibly. If such decays are kinematically impossible or suppressed, then, 
given that couplings between the MSSM and E-sector Higgs fields are generically present 
and imply that the mass eigenstates are mixed, the mostly E-field light Higgses will decay 
into bb, T^T~ or other quarks or leptons depending on the model. Although 15-particles 
would have small direct production cross sections and it would be difficult to detect them 
directly, their presence would be manifest through the dominant Higgs decay modes being 
h 4/, where 4/ symbolically means four SM fields, e.g. bbbb, bbr^r" , r^r^r^r^ , and 
so on ^. This would imply a very complicated Higgs phenomenology. Nevertheless, such 
a scenario is a simple consequence of having an extra sector which couples to the SM or 
MSSM through mixing of the Higgs sectors with one of the extra Higgs mass eigenstates 
being light enough that the SM-like Higgs can decay into a pair them. 

The situation described in the previous paragraph already occurs in the simplest extension 
of the MSSM, the next-to- minimal supersymmetric model (NMSSM) which adds only one 
singlet chiral superfield, S. The very attractive nature of the NMSSM extension of the 
MSSM on general grounds has been discussed for many years j4|; in particular, it avoids the 
need for the /i parameter of the MSSM superpotential term jdHuHd- The NMSSM particle 
content differs from the MSSM by the addition of one CP-even and one CP-odd state in the 
neutral Higgs sector (assuming CP conservation), and one additional neutralino. We will 



follow the conventions of 



51]. Apart from the usual quark and lepton Yukawa couplings, the 



scale invariant superpotential is 



K 



X SH^H, + '-S^ (1) 



^ The situation can be even more complicated if the E-ficld decays into other E-ficlds before the latter 
finally decay to SM fields. In such a case, the SM-like Higgs would effectively decay into 8/. Even more 
complicated variations can also be constructed. 
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depending on two dimensionless couplings A, k beyond the MSSM. [Hatted (unhatted) capi- 
tal letters denote superfields (scalar superfield components).] An effective fi term arises from 
the first term of Eq. ([T]) when the scalar component of S acquires a vacuum expectation 
value, s = (S), yielding 

Hes = As . (2) 
The trilinear soft terms associated with the superpotential terms in Eq. ([1]) are 

XA.SH^Hd + '^A^S' . (3) 

The final input parameter is 

tanP = K/hd, (4) 
where hu = {Hu), hd = {Hd). The vevs hu, hd and s, along with mz, can be viewed as 
determining the three SUSY breaking masses squared for H^, Hd and S (denoted tti^Hu^ 
fn'jj^ and m|) through the three minimization equations of the scalar potential. Thus, as 
compared to the three independent parameters needed in the MSSM context (often chosen 
as /i, tan (3 and M^), the Higgs sector of the NMSSM is described by the six parameters 

\ , K , Ax , A^, tan/5 , . (5) 

(We employ a convention in which all parameters are evaluated at scale rriz unless otherwise 
stated.) We will choose sign conventions for the fields such that A and tan/5 are positive, 
while K, Ax, A,^ and fi^s should be allowed to have either sign. In addition, values must be 
input for the gaugino masses (Mi^2,3) and for the soft terms related to the (third generation) 
squarks and sleptons {rriQ, m\j, m|,, m|, m\, At, At and Ar) that contribute to the radiative 
corrections in the Higgs sector and to the Higgs decay widths. For moderate tan/5, the soft 
parameters which play the most prominent role are m^, mfj, and At. 

Of all the possible new phenomena, the additional Higgses in the NMSSM can lead to, 
perhaps the most intriguing one is the possibility of the lightest CP-even Higgs decaying 
into a pair of the two lightest CP-odd Higgses, hi aiOi, where the latter are mostly 
singlets ^, |9|. Not only would hi aiai decays complicate Higgs searches, but also 

it is found that precisely this scenario can essentially eliminate the fine tuning of EWSB 
in the NMSSM for mn, ~ 100 GeV Q, fl fl- If Br(/ii ^ aiOi) > 0.7 and rua, < 2m^, the 
usual (LEP) limit on the Higgs boson mass does not apply and the SUSY spectrum can be 
arbitrarily light, perhaps just above the experimental bounds, i.e. certainly light enough for 



natural EWSB 
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In this paper, we identify the region of parameter space where this situation occurs 
and discuss how natural this scenario is. The mass of the hghtest CP-odd Higgs boson is 
controlled by the soft-trilinear couplings Ax and and vanishes in the R-symmetry limit, 
Ax, A^ 0. However, both Ax and A^ receive radiative corrections from gaugino masses 
and so arbitrarily small values of trilinear couplings would require cancellation between the 
bare values and the radiative corrections. It turns out that in order to achieve < 2mb 
with Ax, of order the typical radiative corrections, a tuning of trilinear couplings at the 
level of 5% — 10% might be required, but that even less fine-tuning is possible in the context 
of various types of models. Quite interestingly, Br(/ii — > aiOi) > 70 % is automatic in 
this region of parameter space, implying no need for additional tuning in order to achieve 
Br(/ii bb) small enough to escape LEP limits ^. 

It is important to note that to a large extent the tuning in Ax and required for 
Br(/ii — > aitti) > 0.7 and m^^ < can be separated from fine-tuning/naturalness for 
proper EWSB. The value of rriz, as obtained from the renormalization group (RG) equations 
after evolving from the GUT-scale, Mu, down to niz, is primarily sensitive to Ms{Mu), 
m]j^{Mu), ml^{Mu), ml{Mu), ml{Mu), mfj{Mu), ml{Mu), and At{Mu) (and Ab{Mu) if 
tan/? is large). In some cases, M2{Mu) can contribute significantly to the standard measure 
of EWSB fine-tuning with respect to the GUT-scale parameters. In contrast, Mi{Mu), 
Ax{Mu), and Ai^{Mu) have to take very large values in order to contribute significantly to 
fine-tuning. The scenarios that have small fine-tuning that we focus on are ones in which 
these parameters are small enough that they do__not affect the measure of the fine-tuning 
associated with EWSB. In a companion paper 
expanding upon the earlier discussions in Refs. 

Recently various scenarios, similar in spirit, which improve on the naturalness of EWSB 
by modifying Higgs decays or which, in general, suggest the possibility of an extra sector near 
the EW scale have been discussed; see e.g. [ill, [l^lisl. [l^. Phenomenological consequences 



, we discuss EWSB fine-tuning in detail, 



^ Some aspects of EWSB and the possibility of having a hght CP-odd Higgs boson in the NMSSM in the 
R-symmetry hmit were also recently discussed in Ref. [l^l ■ Results and conclusions presented in Ref. [l^] 
are based on studying a region of parameter space very near the R-symmetry limit with soft-trilinear 
couplings much below the typical size of radiative corrections. The conclusions of this reference do not 
necessarily apply to the scenario discussed above that assumes soft-trilinear couplings A\ and of order 
the typical radiative corrections (they could, however, originate from an exact symmetry limit at some 
high scale). We refer to our scenario as NMSSM close to the R-symmetry limit. 
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and possible collider signatures for the NMSSM and, more specifically, for some of these 
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scenarios have also been discussed in Refs. 

The paper is organized as follows. In the next section, we discuss the NMSSM close to the 
R-symmetry limit. Numerical results are presented in Sec. Illll and we conclude in Sec. IIVI 



II. NMSSM CLOSE TO THE R-SYMMETRY LIMIT 

In the NMSSM, one of the CP-odd Higgses is massless in the Peccei-Quinn symmetry 
limit, K — i> 0, or in the R-symmetry limit, A^^, Ax ^ j^. We focus here on the second case 
because it was identified as the easiest way to achieve EWSB without large fine tuning in 
soft SUSY breaking terms j^. This scenario requires that the lightest CP-even Higgs boson, 
hi, decays to a pair of the lightest CP-odd Higgs bosons with BT{hi aiCi) large enough 
that the hi bb signal is highly reduced (compared to the SM signal). Furthermore, it was 
found that this scenario with reduced Br(/ii — > bb) is consistent with the observed excess of 
events at rrih — 100 GeV in Zh production at LEP 7|. 

The masslessness of ai in the limit Afj^,Ax can be understood as a consequence of 
a global f/(l)i? symmetry of the superpotential under which the charge of S is half of the 
charge of HuHd- In the limit A^, — > it is also a symmetry of the scalar potential. This 
symmetry is spontaneously broken by the vevs of Hu-, and S, resulting in a Nambu- 
Goldstone boson in the spectrum. Soft trilinear couplings explicitly break U{1)r and thus 
lift the mass of the ai. For small trilinear couplings, the mass of the lightest CP-odd Higgs 
boson is approximately given as: 



2 ^ I 3AAa cos^ Oa 

ai 



where cos 9a measures the doublet component of the lightest CP-odd Higgs mass eigenstate. 



ai = cos 9a Amssm + sin 9a As- (7) 



In the limit of large tan/5 or |s| ^ v, cos 6*^ can be approximated by 



t> sin 2/3 

cos9a^ -. (8) 



In this limit, the ai mass eigenstate is mostly singlet and 



< - 3. ( '-^^ ^ -^A^] . (9) 



Naively, an arbitrarily small mass for the ai is achievable provided small values of and 
A\ are generated by a SUSY breaking scenario. Indeed, there are SUSY breaking scenarios, 
e.g. gauge mediation or gaugino mediation, which have zero soft trilinear couplings in the 
leading order. However, even if zero values of and Ax are generated at the SUSY breaking 
scale, the corresponding electroweak (EW) scale values will be shifted due to radiative 
corrections from gaugino masses. The typical EW scale values can be estimated from the 
one- loop renormalization group equations for A\ and A^, 



dA 



X 



dt 
dA, 



167r2 
12 



6AAt + SX^Ax + 4k^A, + 6glM2 + (^^ ] gfM, 



(10) 

dt 167r2 ^" ' ^^^^ ' ^^^^ 
where t = \og{Q/mz), At is the top soft trilinear coupling and M2 and Mi are the masses of 
the SU(2) and U(l) gauginos, and we have neglected terms proportional to Ab and Ar. Start- 
ing with A\ = at the GUT scale, we find that A\{mz) ~ M2 because the loglMu/mz) 
coming from the integration approximately cancels the (6(7|)/(167r^) loop factor. On the 
other hand, receives contributions from gaugino masses only at the two-loop level imply- 
ing that A^{mz) is expected to be much smaller than Axlrriz). Assuming gaugino masses of 
order 100 GeV, we should naturally expect A\{mz) — 100 GeV and Anijnz) ~ few GeV. 
Much smaller values would require cancellations between the values of Aa, coming from 
a particular SUSY breaking scenario and the contributions from the radiative corrections. 

For sizable A\{mz), Eq. ([9]) is no longer a good approximation for the mass of the 
lightest CP-odd Higgs. In order to understand the numerical results presented later, we 



have developed a more accurate formula. The 2x2 mass matrix squared 
Higgs bosons, in the basis {Amssm, As), has the following matrix elements 



Ml 
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Ml, 



The eigenstate masses are 

ti2 Z 



2As 
sin 2/3 
{Ax 



KS] 



2ks) 



2\kv^ sm2p + \A, 



f ^ sin 2/9 
^ 2s 



SkA/^s . 



{Ml, - Mi,y + 4{Ml,y 



The mixing angle for the diagonalization process is obtained from 

sin26'A = 



2Mf, 



^{Ml-Ml^^ + A{Ml^^ 



'or the CP-odd 

I 

(12) 
(13) 
(14) 

(15) 

(16) 



where we are using the convention defined in Eq. ([7j). Obviously, the value of ^A_is only 
determined mod(-7r). In our numerical work, we employ the program NMHDECAYjSl] which 
adopts the convention < 6*^ < vr. If Mf^ is much larger in magnitude than the other 
entries, then it must be positive and the mass squared of the lightest CP-odd Higgs boson 
is then given by 

For typical values of trilinear couplings, \Af^\ <^ \Ax\ ~ f ^ and tan/3 > few we find: 

ml ~ 3s , 19) 

which reduces to Eq. ([9]) if we neglect A\ and compared to s. Similarly, the mixing angle 
is determined by 

cob29, = 2cob'^9,-1^I^|M|__i. (20) 

sin 29a = 2 sin 9a cos 9a ^ -jt^-^jjt > (21) 

from which we obtain (using the conventions defined earlier) sin6'A ~ 1 and the doublet 
component of Oi is given by 



eo,g,^ ^ HA,-2ns)sin2(5 

Mf^ - MI2 2Xs{Ax + Ks) + 3kA^s sin 2/3' ^ ^ 



This reduces to Eq. ([8]) for very small Ax and A^. 

Both Eq. and Eq. ( |T9l) indicate that sizable Ax could give a large contribution to the 
mass of the lightest CP-odd Higgs. However this term is highly suppressed if the lightest CP- 
odd Higgs is mostly singlet. In this case, both terms in Eq. (Q or Eq. (fT9|) are comparable 
and then it depends on their relative sign as to whether they contribute constructively 
or destructively. One measure of the tuning in Ax{mz) and Af^{mz) necessary to achieve 
ma^ < 2mb is 

FMAx = max{\FA,l\FAj} , (23) 

where 

_ Axjmz) dml^ _ AJjnz)dm^ 



are evaluated for given choices of Axlrnz) and A^imz) that yield a given . This definition, 
which reflects the fact that rria^ is determined by these parameters for fixed A, k, /icfr, tan/3, 
will prove useful below for discussing the sensitivity of nia^ to GUT-scale parameters. We 
argue below that Fmax is typically an upper bound on the magnitude of fine tuning with re- 
spect to GUT scale parameters. Thus, Fmax provides a useful first measure for determining 
a "preferred" region of parameter space where a small mass for the lightest CP-odd Higgs 
boson is achieved with the least tension. However, we will also find that in large classes of 
models Fmax greatly over-estimates the fine-tuning. 

The fine-tuning measure for m^^ relative to GUT scale parameters that is completely 
analogous to that employed for EWSB is 

d log Tfl^ 

Fma, = "^^^pfp , with fp = — ^ , (25) 

where p is any GUT-scale parameter, 

p = Mi{i = 1,2,3), mjj^, mjj^, m|. At, A^, Ax, rrig, mfj, mjj (26) 

(to name the most important GUT-scale parameters). However, for our purposes we can 
simplify this general computation because we are only interested in m^^ fine tuning for cases 
in which the EWSB fine tuning is already known to be small, that is cases in which 



rrir 



-A {mz)s imz) + 



(27) 



tan^ P{mz) — 1 

is insensitive to the parameters p. Small fine tuning for m| means that v{mz) (which sets 
rriz), s{mz), tan Plruz), ni\j^{mz) and m?^^{mz) are not fine-tuned with respect to the 
various p listed above. The only additional parameters upon which m^^ depends that could 
still be sensitive to the GUT-scale parameters p when niz is not are A\{mz) and A^^mz). 
For any of the p, we can then approximate 

/ P F dAxjmz) , p ^ dA^jmz) . . 

Ax{mz) dp A^{mz) dp 

To understand the implications of this formula, we need to solve the RG equations and 
express Ax{mz) and Ai^{mz) in terms of the GUT-scale values of all the soft-SUSY-breaking 
parameters. (Of course, we are taking the GUT scale as an example; a similar exercise can 
be done for any scale.) The solution depends on A, k and tan /5; we give only a representative 
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example. For A = 0.2, k = ±0.2 and tan/9 = 10 we find: 



Aximz) = -0.03A^ + 0.93Aa - 0.35A - O.OSAfi - 0.37M2 + O.66M3, (29) 
A^{mz) = 0.90A«-0.11Aa + 0.02A + 0.003Mi + 0.025M2-0.017M3, (30) 

where the parameters on the right-hand side of these equations are the GUT-scale values. 

Before discussing the fine-tuning implications of Eqs. fl29l) and (l30l) . it is useful to un- 
derstand a few of their features. First, we note that while the coefficients in the Ax{mz) 
expansion do not change much when changing A and n, the coefficients in the A^{mz) ex- 
pansion are quite sensitive to changes of the A, n Yukawa couplings. The reason is that 
gaugino masses enter through the Ax term in the RG equation, see Eqs. (fTU]l - (fTTl) . the 
strength of which is controlled by A. The opposite sign of the coefficient in front of M3 in 
the expansion of Ax as compared to the coefficients in front of Mi and M2 is due to the 
fact that M3 enters only through A^ in the RG evolution, while both Mi and M2 enter 
directly. Similarly, the opposite sign in front of the gaugino masses in the Ax expansion as 
compared to the expansion is due to the fact that gaugino masses enter A^ through Ax 
in the evolution. The above expansions provide several ways to achieve A^ Ax ^ mz- 
The easiest way is to assume that a SUSY breaking scenario generates negligible trilinear 
couplings at the GUT scale. 

We now return to a consideration of the fine-tuning implications of Eqs. (129|) and (!30ll . 
If there is a px that dominates Ax{mz) and a that dominantes A^{mz), and these are 
different, then 

^^^^"^^^-(9(1), and -r^^^4^-0{l), (31) 



Ax{mz) dpx ' A^,{mz) dp^ 

and roughly F^^^ ~ Fmax- If the same p dominates both Axijnz) and A^i^rriz) then 
Fma-^ ^ fp ^ Fa, + Fa^- This result also holds if the GUT-scale parameters are correlated. 
For example, consider the case of universal gaugino masses and zero trilinear couplings at 
the GUT scale, for which Ax{mz) = O.26M1/2 and A^{mz) = O.OIM1/2. Then, 

Ml/2 dAx{mz) Ml/2 dA^{mz) 



Ax{mz) rfMi/2 A^{mz) rfMi/2 
and it is quite precisely the case that 



(32) 



Fm.^ = fM,,, ~ Fa, + Fa^ (33) 
10 



As we shall see in the numerical section of this paper, Fa^ and Fa^ are typically opposite in 
sign and of similar magnitude. This can already be seen from the approximate formula of 
Eq. ^ where the linearity of m^^ in Ax{mz) and A^{mz) would give (neglecting the mild 
dependence of s{mz), ... on M1/2 for the points of interest) 

FA, + FA^r^l, (34) 

and fine tuning with respect to M1/2 would be small. The more precise (but still approx- 
imate) result of Eq. flTI?]) gives somewhat larger results for Fa^ + Fa^ for some parameter 
choices, but a subset of choices still gives Fa^ + Fa^ ~ (!)• Later, we will present numer- 
ical results for Fa-^ and Fa^ that confirm that they largely cancel against one another for a 
significant fraction of parameter choices. The result of Eq. fl33l) clearly applies whenever the 
gaugino masses are correlated (in any way) and trilinears (at the GUT scale) are small or 
correlated with the gaugino masses. As already noted, this same result also applies whenever 
a single term dominates in Eqs. (12^ and (150]) . Many models fall into one or the other of 
these categories. 

Let us reemphasize that tuning in A\{mz) and Ai^{mz) is completely unnecessary to 
achieve a light CP-odd Higgs boson in models with specific relations among GUT-scale 
parameters. Any SUSY breaking scenario that determines all soft trilinear couplings and 
gaugino masses from a SUSY breaking scale will automatically give A\{mz) = c\Msusy and 
A^i'mz) = c^MsusY, where cx and depend (given the known values of gi and (72) only on 
the couplings A, k, and (equivalently, tan/3, given the known value of my/). The mass of 
the lightest CP odd Higgs boson will be given as m^^ = /(A, k, tan /3)M|^5y and will either 
be small or not. This means that, in any SUSY breaking scenario that is determined by 
a SUSY breaking scale only, there is no tuning of with respect to the SUSY breaking 
scale. (Algebraically, ,,'^'°f,T"^ = L) This result holds even if there are large cancellations 
among the RG contributions to A\{mz) and Af^{mz)- Whether or not a light CP odd Higgs 
boson is possible simply depends on the above couplings. 

Let us now turn our attention to 'Bi{hi cticti), a completely general expression for 
which is the following (neglecting phase space suppression): 



+ 0^ 

— COS 2(3 sin(/3 -|- a) cos^ 6a cos 63 



327rg^mh^ 

+2A^ ( sin(/3 — a) sin^ 6a cos 63 -\ — cos^ 6a sin 6 
V V 
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+ (cos a siv? /3 — sin a cos^ [3) cos^ 6a cos 6^ 

+2\k ^cos(/5 + a) sin^ 6a cos 6*5' — sin 26a sin 6s 

s s 
— sin(/5 — a) sin 26a cos 6'^ H — sin 2(3 cos^ 6*^ sin 6*5 

V V 

+— — - (sm(/3 — a) sin 26a cos 6s + sin 2/5 cos^ 6a sin 6'5) 

/4k2s 2/€yl,\ .2 ■ " 

+ sm 6a sm 6s 

V V I 



(35) 



This is to be contrasted with the dominant SM decay channel, hi — > 66, the width for which 
is (neglecting phase space suppression): 

T{hi^hh) ~ ^TT^^ ( — ^ ) rrih.ml cos^ 6s , (36) 



327rm^ \cos/3^ 

Our conventions are those of NMHDECAY in which the WW coupling relative to SM 
strength is given by sin(/3 — a) cos 6*5 and v"^ = 2m^/(7|; is to be evaulated at scale rrih^. 
In Eq. (I7j), cos6'a gives the MSSM doublet component of the ai and sin6's is the coefficient of 
the singlet component of the hi, both are small numbers in our scenario. Detailed numerical 
results and discussion for the hi — > aiai width will be given in the next section. However, 
it is important to understand the limit in which A\,Ak 0. This is nicely illustrated in 
the case of v/s <^ 1. In this limit we have: a — ^ /5 — 7r/2 (more precisely, sin(/5 — a) — *• 
sm2(3 + O (vys^)); cos6a ^ (v/s) sm2(3 + O (vys^) (see Eq. (IHD); sin^A ^ 1 - C (vys^); 
sm6s -> -{v/s){Xy{2K^)){l-{K/X)sm2(3) + 0{vys^); cos6s l-0{vys^). With these 
inputs, we find 

. mlr r^,o ^, . ^„ 4:Xks „ 4k^s 

i [hi — >• aiOij 



32'iTg^mh, 



m 



w 



2\^ + 2\k sin 2/5 cos 6a + sin 6s + iv"^ I s^) 

V V 

[0{vVs')Y. (37) 



327rg^mh^ 

In particular, the [. . .] does not actually vanish for Ax = = (there is no exact symmetry 
argument). One finds that Br(/2,i aiai) can approach the [0.1 —0.2] range for the smallest 
K values allowed by general theoretical consistency (good EWSB vacuum, . . . ) for a given 
A value. This is, of course, insufficient for escaping LEP constraints when ruhj^ ~ 100 GeV. 
Thus, having an adequate size for Br(/2,i — > aiOi) for escaping the LEP constraints depends 
critically upon having non-zero values for A^ and, in particular, A^. From the numerical 
results presented in the next section, it will be clear that magnitudes for Ax and A^ of order 
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those developed from RGE running beginning with Ax{Mij) = Ai^{Mu) = are sufficient to 
give large Br(/ii aiOi). 

III. RESULTS 

As discussed in the introduction, the Higgs sector in the NMSSM is determined by 6 
basic parameters, given in Eq. (jl]), as well as subsidiary parameters entering through loop 
corrections. (In this section, all parameters are defined at scale mz-) Consequently, a 
complete survey of the parameter space is difficult. To present results in a manageable way, 
we fix /i and tan /5 together with all soft SUSY breaking masses and scan over trilinear and 
soft-trilinear couplings. We will plot results in various 2-parameter planes. The parameters 
are scanned over the following regions with fixed steps: A G (0, 0.5) using 30 steps of size 

0. 01666; K E (-0.5, 0.5) using 70 steps of size 0.014286; Ax E (-300 GeV, 300 GeV) using 100 
steps of size 6 GeV; and finally A^ E (-20 GeV, 20 GeV) using 100 steps of size 0.4 GeV. 
Varying the fixed soft SUSY breaking masses does not have any significant effect on the 
results while different choices of /i and tan/3 lead to important changes. Thus, we will 
present results for several choices of /i and tan /5 keeping all SUSY-breaking masses fixed at 
MsusY = 300 GeV. Let us recall that we are interested in looking for parameters yielding 
m/j^ > 90 GeV and m^^ < 2m;,. The latter is imposed so that the LEP limits on hi aiai — »• 
46 do not apply (since if the ai decayed primarily to bb they would require rrih^ ^110 GeV, 

1. e. above our preferred rrih-^ ~ 100 GeV value). The m/^^ > 90 GeV restriction implies that 
we are above the maximum value for which LEP limits on hi aiOi At are available. 
In our plots, the small dark blue diamonds are all points that satisfy the above constraints, 
while the large light blue crosses are those which satisfy all experimental limits, the main 
experimental constraint being that Br(/ii bb) must be suppressed sufficiently by a large 
Br(/ii — s> oiOi) that LEP limits on the Zhi Zbb channel are satisfied. Roughly, this 
requires Bi{hi aiOi) > 0.7. 

We first present results for Msusy = 300 GeV, ^ = 150 GeV and tan/5 = 10. In Fig. [1] 
we plot the allowed region of parameter space in the A^ — Ax and k — A planes. Similarly, 
in Fig. [2] we plot the allowed region in the A^ — k, and Ax — n planes. In Fig. [HI we plot 
a selection of the range of and Ax values that have < 2m^ for fixed values of A 
and K. From Fig. [3l we see that for a given value of Ax{mz) keeping m^^ < 2mb requires 
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that A^imz) be adjusted to a level of order 10%; at fixed A,^{mz)i A\{mz) must lie within 
about a 5% range. As quantified later, this is a rough measure of the tuning required to 
achieve consistency of the envisioned scenario with LEP constraints when m^^ ~ 100 GeV. 
Note that A^, Ax tuning could be made arbitrarily mild if very small values of Af^{mz) were 
allowed. However, as illustrated by the dark blue points, too small a value for A^{mz) leads 
to a value for Br(/ii — s> aiOi) that is small and therefore a value of Br(/ii — bb)) that is too 
large for consistency with LEP constraints on the Zh Zbb channel. In any case, as we 
have discussed, very small values of Af^{mz) and Ax{mz) would be purely accidental from 
the RGE point of view. 

The correlations between Ai^{mz) and Ax{mz) required to achieve m^^ < can be 
understood from Eqs. ffT^ - f|T9]) . For the following discussion, we consider the case of /i > 0, 
implying s > 0. Clearly, in order for both eigenvalues of the CP-odd Higgs mass matrix 
squared to be positive, both diagonal elements have to be positive. In the limit \AJ^ <^ 
\Ax\ ^ V <^ s, and tan/3 > few, is dominated by the last term (unless \A^^^ is too 
small), and so Mfg > leads to the condition < 0. This means that the contribution 
of the nAt^ term to the mass of the lightest CP-odd Higgs is positive, see Eqs. ([9]) and ( |T9l) . 
Next, we note that Mf^ > requires Ax + us > 0. Given that nA,^ < 0, this implies that 
the denominator of the first term in the m^^ expression given in Eq. (fT9|) is positive, and 
that this first term will therefore have the same sign as nAx- Thus, in order for the terms 
proportional to kAx and —nAi^ in Eq. ( fT9l) to cancel so as to give small ?Ti^^, kA\ must 
have the same sign as nA,^. Given that < 0, we must therefore also have nAx < 0. 
Altogether, we have two disconnected regions of allowed parameter space. The first one is 
for K > 0, < and < 0. It is the largest region because the condition -|- > is 
easily satisfied, especially with smaller values of Ax. The second region is for k < 0, > 
and > which is further constrained by -|- > and thus requires larger values of 
Ax and consequently larger values of A^. 

The above discussion is not valid when is so small that the kA^ term does not dominate 
MI2. For such parameters, < 2mf, can still be achieved, but Br(/ii aiai) > 0.7 is 
not possible. This region is seen in Fig. [2] as a large dark (blue) region with a large range 
of Ax > with K > in the right-hand plot and a narrow band of very small A^ values 
and K > in the left-hand plot. For this region, A < 0.1 is typical. This region can be 
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Mj^=300 GeV, /t=150 GeV. tan|S=10. mi_>90 GeV. m^^KZia^ 
allowed by LEP 




Mbiiby=300 GeV, /i=150 GeV, tani8=10. m^^>90 GeV, m„<2mj 

X = allowed by LEP 



0.4 
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0.0 



-0.2 
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FIG. 1: Allowed parameter space in the A,^ — A\ and k — X planes. Light grey (cyan) large crosses 
are points that satisfy all experimental limits. The dark (blue) diamonds are those points that do 
not have large enough Br(/ii — > aiai) so as to suppress Br(/ii — > 66) sufficiently to escape LEP 
limits. 



understood by noting that when kA^ is negligible, then Eq. (fT9!) reduces to 

2 A^-^o OAxkXv"^ sm2P 



m 



(38) 



2{Ax + Ks) 

Since we require both m^^ > and Mf^ oc {Ax + ks) > 0, see Eq. ( IT2|) . we must have Ax > 
and K > 0. Small m^^ is easily achieved in a variety of ways, for example if Ax is small 
compared to ks or, more generally, for small XAx and large tan/?. Radiative corrections 
play a significant role also, reducing the tree-level prediction given above by a substantial 
amount. 

The dependence of the branching ratio for hi — ^ aiOi on and Ax is given in Fig. HI As 
we have stressed, a large value for Br(/ii — > aiOi) is crucial for an rrih^ ~ 100 GeV light SM- 
like Higgs to have escaped LEP constraints. (It is the constraint m^-^ < 2mb which guarantees 
that when ~ then XAx must be very small as well, and vice versa, thereby implying 
small Br(/ii aiOi).) It is interesting to see that if \Af^\ > 2 GeV and \Ax\ > 30 GeV (as 
typical for the sizes of the RGE-induced contributions) then Br(/ii — > aiai) is almost always 
large enough for the hi to have escaped detection through the hi —>■ bb channel. These plots 
also show clearly that Br(/ii aiOi) approaches a small (most typically, extremely small) 
value when Ai^,Ax 0. This suppression was discussed analytically in the small v/s limit 
in the previous section. 
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MBtiBr=300 GeV, /i=150 GeV. tan(S=10, m^>BO GeV, m„<2mj MBtrar=300 GeV, /i=150 GeV, tan(S=10, inij>90 GeV, m„<2mj 
X = allo-wed by LEP x = allowed by LEP 




FIG. 2: Allowed parameter space in the — k and Ax — k, planes. Point conventions as in Fig. [H 



MsnsY=300 GeV. At=150 GeV, tan/S=10. X=0.38, K=0.4, m^>90 GeV 




A^(mz) (GeV) 



FIG. 3: A selected region of the allowed parameter space in the — A\ plane for fixed values of 
A = 0.38 and k = 0.4. Point conventions as in Fig. [TJ 

Another way of characterizing the hmited allowed region in and Ax for which 
stays below 2mb (shown in Fig. [3]), is to calculate the tuning [as defined in Eq. ( 1231) ] of Ax 
and Ai^ necessary to achieve this situation. The dependence of Fmax on and Ax is given 
in Fig. [51 Similarly, the dependence of Fmax on k and A is given in Fig. [61 As expected, the 
smallest fine tuning or sensitivity is achieved for as small A^ and Ax as possible. Of course, 
as we discussed earlier, very small values of and Ax would require cancellations between 
the bare values and the RGE-induced radiative corrections and this kind of cancellation 
would not be visible from the definition of Fmax given in Eq. ( l23l) . However, this is not a 
particularly worrisome point given that, as discussed with regard to Fig. [H very small values 
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Mj^=300 GeV, /t=150 GeV. tan|S=10. Tn^>90 GeV. m,^<Zm^ 



allowed by LEP 




A. (GeV) 



FIG. 4: Br(/ii — >■ aiai) vs. and A\ for ^ = 

Fig.m 

Mj^=300 GeV, /i=150 GeV. tanjS=10. mi_>90 GeV. m,^<Zm^ 




A. (GeV) 



Mj^=300 GeV, /t=150 GeV. tan|S=10. Tn^>90 GeV. m,_<2in5 
allowed by LEP 




-200 200 

Ax (GeV) 

150 GeV and tan (3 = 10. Point conventions as in 



Mj^=3O0 GeV, /i=150 GeV. tanjS=10. mi_>90 GeV. m^^KZia^ 




-200 

Ax (GeV) 



FIG. 5: Tuning in m^^ vs. and Ax- Point conventions as in Fig. [H 

of Ai^ and Ax do not in any case lead to large enough Br(/ii aiOi) that Br(/ii — > bb) is 
adequately suppressed. Thus, in the region of parameter space where soft trilinear couplings 
are at least of order the typical RGE-induced contributions, which is also the region where 
Br(/?,i ttitti) is large, the Fmax measure of the tuning of the soft trilinear couplings can 
be as small as C(5% - 10%). 

However, as discussed in the analytic section of the paper, fine-tuning with respect to 
GUT-scale parameters, denoted there as Fm^^, need not be as large as Fmax- In particular, 
there are many model scenarios in which Fm^-^ is reduced compared to Fmax by cancellations 
between the dependence of Ax{mz) on the dominant GUT-scale parameter p and the de- 



17 



100 
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X = allowed bv I,KP 
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FIG. 6: Tuning in vs. k and A. Point conventions as in Fig. [H 
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FIG. 7: In the left-hand frame, we plot F^^ vs. F^^ for the points with Fmax < 100. In the 
right-hand frame, we plot F^^ + Fa^ vs. A^^ for these same points. 



pendence of Ai^{mz) on this same p, in the simplest cases yielding Fm^^ ~ Fa^ + Fa^- In the 
most naive approximation, m^^ is linear in Ax and linear in and therefore Fa^ ~ ~Fai^ 
and -Fma;^ is then quite small. Thus, it is important to understand the correlations between 
and Fa^ in order to assess the extent to which -|- Fa^ can be small. In the left-hand 
window of Fig. [3, we plot Fa^ vs. F^i^ for points with Fmax < 100. We see that many of the 
points have F^i^ ~ —Fa^- The corresponding values of -|- Fa^ appear in the right-hand 
plot. We observe that the LEP-allowed points with small negative (which are those 
with the smallest magnitudes of and Fa^ separately) yield still smaller Fa^^ + Fa^ and 
therefore possibly very small F^^^ for appropriate GUT-scale models. For all LEP-allowed 
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FIG. 8: In the left-hand frame, we plot Fmax vs. for the points with Fmax < 100. In the 
right-hand frame, we plot F^^ + F^^ vs. for these same points. 



Mj^=300 GeV, /i=150 GeV. tanjS=10. mi_>90 GeV. m^^KZia^ 

X = allowed by LEP 



Mj^=3O0 GeV, /i=150 GeV. tanjS=10. mi_>90 GeV. m,^<2m^ 

X = allowed by LEP 




FIG. 9: Singlet component of hi vs. and Ax. Point conventions as in Fig. [H 
< points, Fa^ + Fa^ is much smaller than Fmax- 

It is also useful to understand the extent to which Fmax and Fa^ + Fa,, depend on nta^ . 
This dependence is revealed in Fig. |H1 We see that the smallest Fmax values are achieved for 
the larger values up near 2mfe. Small values of Fa^ are distributed over the broader 
range of 4 GeV < nia^ < 2m;,. In either case, the fine-tuning associated with m^^ < 2mb 
suggests some preference for m^^ > 2m^. This is important in that hi ^ aiOi 4r decays 
may prove to be visible at the LHC, whereas hi aiOi with ai cc,gg, . . . will be much 
harder to detect. 

Higgs phenomenology is crucially dependent upon the singlet and doublet compositions of 
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Mj^=300 GeV, /i=150 GeV. tanjS=10. ini_>90 GeV. m,^<Zm^ Mj^=3O0 GeV, /i=150 GeV. tanjS=10. ini_>90 GeV. m,^<Zm^ 

allowed by LEP x = allowed by LEP 
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FIG. 10: Singlet component of ai vs. and Ax. Point conventions as in Fig. [H 



Mj^=300 GeV, /i=150 GeV. tanjS=10. mi_>90 GeV. m^^KZia^ 




&! singlet component 



FIG. 11: Singlet component of hi vs. singlet component of ai. Point conventions as in Fig. [H 

the lightest CP-even and CP-odd Higgs mass eigenstates. The hghtest CP-even Higgs boson 
is preferably SM-like although it can have up to 20 % singlet component, see Fig. [H while 
the lightest CP-odd Higgs boson has to be very close to being a singlet, at least 98 % singlet 
in the region where the least tuning is necessary, see Fig. [lOl (Here, sin6'yi, see Eq. ((71), 
is the singlet component of the ai at the amplitude level. The probability for the ai to be 
singlet is sin^ 9a- Similar definitions are used in the case of the hi.) This correlation between 
the hi and ai compositions is explicit in the plot of the hi singlet component versus the ai 
singlet component given in Fig. [TT] This figure shows that as the Oi becomes less singlet 
the hi must have less singlet component. The left-hand plot of Fig. [12] also makes it clear 
that the hi singlet component must be small (or zero) in order to maximize Br(/ii — > aiOi). 
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Mj^=300 GeV, /t=150 GeV. tan|S=10. m^>90 GeV. m^^KZia^ 
allowed by LEP 
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allowed by LEP 
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FIG. 12: Br(/ii — > oiai) vs. singlet component of hi and ai. Point conventions as in Fig. [TJ 
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FIG. 13: Br(/ii aiai) and '^s. the non-singlet component of the ai. Point conventions 

as in Fig. [TJ 



The right-hand plot of Fig. [12] shows that the ai singlet component must be at least 70% 
in amplitude, i.e. 50% in probability, for the range of parameters scanned with most points 
corresponding to an Oi that is mainly singlet. 

A particularly revealing pair of plots are those of Fig. [13] where we show how Br(/ii 
aifli) and Fa^ + Fa^ depend on cos6'a, the non-singlet component of the ai. We see from 
the first plot that there is a lower bound on | cos 6*^1 of order 0.06 (this is the bound for 
tan/3 =10 — the precise number depends on tan/5) in order for 'Bi{hi aiOi) to be large 
enough that LEP limits are evaded. The existence of the lower bound follows from the fact, 
discussed earlier, that k and A\ must have opposite signs in order for the scenario to be 
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FIG. 14: Allowed parameter space in the A^ — A\ and k — A planes for tan/3 = 3 and fi = 150 GeV. 
Point conventions as in Fig. [TJ 

viable. Since, the non-singlet component of the ai is proportional to Ax — 2ks, see Eq. (l22l) . 
it cannot go to zero for the given correlation of signs. From the second plot, we see that the 
Fa^ + -Faa measure of fine tuning for m^^ has a very distinct minimum value close to for 
cos 9a ~ —0.1. Combined with the dependence of Fa^ + Fa^ on displayed in Fig. [8], we 
see a possible preference for > 2m^ and cos6'a ~ —0.1 in order to be reasonably certain 
that fine-tuning is not required in order to achieve large Br(/ii — > aiOi) and m^^ < 2mh. We 
note that the coupling of the ai to down type quarks is proportional to tan (3 cos 6a which 
is never smaller than about 0.6 in magnitude for the points that escape LEP limits. Thus, 
even though the ai is largely singlet, it has substantial down-quark and lepton couplings. 
In a follow-up paper 



19|, we show that this implies an ma^-dependent lower limit on the 
branching ratio for T — 701 decays. This lower limit can probably be probed at future, 
if not existing, B factories if is not too close to Mx, especially if m^^ > 2771^- so that 
fli T^T" decays are dominant. 

Finally, in order to see the effect of varying tan (3 and /i, we give a selection of some of the 
same plots for: /i = 150 GeV and tan /5 = 3; /i = 150 GeV and tan /3 = 50; and /i = 400 GeV 
with tan/5 = 10. These appear in Figs. [T3] — [23l We note that there is an exact symmetry 
under /i —fi (implying s ^ — s in the A > convention we employ), —A,^ and 

Ax — > —Ax- Thus, only positive fi values need be considered when both signs of Ax and 
are scanned. 

Comparing the plots for tan/? = 3, Figs. [H] — [16], and the plots for tan/? = 50, Figs. [17] 
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MBtiBr=300 GeV, /i=150 GeV. tan(S=3, mij>90 GeV. m^^<2m^ 

X = allo-wed bv LEP 
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FIG. 15: Allowed parameter space in the — k and ^4^ — k planes for tan/3 = 3 and ^ = 150 GeV. 
Point conventions as in Fig. [TJ 



with the previous plots for tan/3 = 10 we see that the region of allowed parameter 
space for k > 0, < and < does not change much when varying tan/3. On the 
other hand the allowed region of parameter space with k < 0, > and > changes 
dramatically. It does not even exist for tan/5 = 3 ^. The reason is that the first two terms 
in Eq. f|T^ are less suppressed for smaller tan (3 and they combine into a term proportional 
to Ax + 4ks. The last term in Eq. dominates only for larger values of A^. If the last 
term does not dominate, the first two terms are positive only for very large Ax making the 
region with k < 0, > and > very constrained. For exactly the same reason, this 
region is larger for tan/3 = 50 compared to tan/3 = 10 and the allowed parameter space 
expands to lower values of as particularly noticeable in the vs. k plots; compare 
Figs. [18], [T5l and [2l In the region of parameter space with k > 0, < and < the 
condition Ax + > is always satisfied when the necessary condition for having positive 
mass squared eigenvalues. Ax + ns > 0, is satisfied. Therefore, this region of parameter 
space is not very sensitive to tan /3. 

Quite interestingly, once < 2mfe is achieved for tan/3 = 3, Br(/ii — ^ aiOi) is au- 
tomatically large enough (> 0.75 for the most part) and, correspondingly, Br(/ii — * bb) is 
sufficiently suppressed that Zh Zbb constraints at LEP are satisfied; see Fig. [161 As in 



^ There is a tiny allowed region with positive A^, and positive but very small A\. This region is not a 
continuation of the region present in the plots for tan P = 10. 
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FIG. 16: Fmax and Br(/ii — > aiai) vs. A,^ for tan/3 = 3 and /i = 150 GeV. Point conventions as 
in Fig. [H 
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FIG. 17: Allowed parameter space in the Af^ — Ax and k — X planes for tan/? = 50 and /x = 150 GeV. 
Point conventions as in Fig. [TJ 

the tan (3 = 10 case, Fmax is smallest for the smallest possible and approaches ~ 15 for 
a couple of the points with > 0. Overall, it would appear that rriaj^ < 2mh is a bit more 
difficult to achieve without tuning when tan/5 = 3 vs. when tan/5 = 10. For tan/5 = 50, 
see Fig. [191 small Fmax 20) is achieved for a narrower range of and > gives 
the smallest Fmax values (whereas for tan/5 = 10 there was no particular preference with 
regard to the sign of A^). 

Finally, we discuss the effect of increasing fi. We consider /i = 400 GeV at tan /5 = 10 and 
present results in Figs. [20] — [231 These figures clearly show that the range of parameter space 
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FIG. 18: Allowed parameter space in the A,^ — k and Ax — k planes for tan /? = 50 and fi = 150 GeV. 



Point conventions as in Fig. [TJ 
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FIG. 19: Fmax and Br(/ii aiai) vs. for tan/3 = 50 and = 150 GeV. Point conventions as 
in Fig. [TJ 



for which m^^ < 2mfc shrinks with increasing /i. This is easy to understand from Eq. ([9]) 
or Eq. f|T9l) . For fixed A, increasing /i results in an increase of s. Consequently the term 
proportional to A\ in the formula for m^^ is further suppressed while the term proportional 
to is enhanced. In order to compensate for this so as to keep small, smaller values of 
K and A^ and larger values of A and A\ are required as compared to the /i = 150 GeV case. 
Moreover, the < 2mb region of parameter space with k < 0, > and > is even 
further constrained by the condition A\ + his > which requires larger A\ for larger s. These 
effects are clearly visible in Fig. [20] and Fig. [2TJ The increased value of Ax required for this 
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Mj^=300 GeV, /t=400 GeV. tan|S=10. mi_>90 GeV. m^^KZia^ 
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FIG. 20: Allowed parameter space in the Ai^ — Ax and k — X planes for tan /? = 10 and fi = 400 GeV. 
Point conventions as in Fig. [TJ 
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FIG. 21: Allowed parameter space in the Af^ — K and Ax — k planes for tan/3 = 10 and /x = 400 GeV. 
Point conventions as in Fig. [TJ 

scenario to work for larger fi leads to larger tuning necessary to achieve < 2mb. Fig. [22] 
shows a dramatic narrowing of the A^^ region for which moderate Fmax values are achieved 
and that the best Fmax value for points having large enough Br(/ii aiOi) that LEP limits 
are evaded is of order ~ 30 in contrast to the Fmax ~ 15 values achieved for /i = 150 GeV. 
However, we see from Fig. [23] that the values of and Fa^ are more strongly correlated 
so that Fa)^ + Fa^ can again take on small values. As we have seen, the latter implies small 
Fma^ in models in which the soft-SUSY breaking parameters are correlated or A\{mz) and 
Af^iniz) are dominated by a single soft-SUSY-breaking parameter. 
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FIG. 22: Fmax and Br(/ii — > aioi) vs. for tan/3 = 10 and /i = 400 GeV. Point conventions as 
in Fig. [H 
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FIG. 23: We plot Fa, vs. Fa^ (left window) and Fa„+Fa^ vs. (right window) for /i = 400 GeV, 
and tan (3 = 10. Point conventions as in Fig. [TJ 

Given that detection of the hi at hadron colliders will be quite challenging for the large- 
Br(/ii — > aiCi) scenarios that we focus on, an interesting question is then whether or not the 
other Higgs bosons /i2,3 and 02 might be detectable. This will depend on their masses and 
on their couplings. We will see that within the scenarios considered, their masses can range 
from somewhat above 100 GeV to quite large values. As regards their couplings, since the 
hi has quite SM-like WW, ZZ couplings for the scenarios being considered (as illustrated 
earlier), the /;.2,3 will have quite weak WW, ZZ couplings, and of course the a2 has no tree- 
level couplings to WW, ZZ. Thus, the question is whether production mechanisms relying 
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FIG. 24: In the left-hand set of plots, we present rrih^ vs. rn-aj (top), m/^,, vs. (middle) and 

m/i3 vs. m/j2 (bottom). In the right-hand set of plots, we present Cd{a2) vs. rua^ (top), Cd{h2) vs. 
m/i2 (middle) and CdQiz) vs. (bottom). All plots are for tan/3 = 50 and ^ = 150 GeV. Point 
conventions as in Fig. [TJ 
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FIG. 25: We plot Cdih^) vs. Crf(/i2) for tan/3 = 50 and = 150 GeV. Point conventions as in 

Fig.m 

on qq couplings of the /i2, ^3, 02 could lead to an observable signal. This same question of 
course applies to the H and A of the MSSM. The answer there is that high tan (3 is required, 
in which case gg bbH, bbA production is highly enhanced since the H, A have bb coupling 
strength of order tan f3 times the SM-like strength. We wish to explore the extent to which 
this also applies in the NMSSM for the h2, h^, a2 in scenarios with large Br(/ii aiai). For 
this purpose, we present some additional plots in the case of tan f3 = 50. We will denote the 
strength of the bb coupling (generally any down-type quark or lepton) of any given Higgs 
boson relative to the SM-like strength by Cd{h), where h = a2,h2,h3 will be considered. 
The relevant plots appear in Fig. [2H 

In the left-hand plots of Fig. [211 we give nih^ vs. rria^ , nih^ vs. and m^g vs. rrih^ ■ The 
light (cyan) points are those that have both m^j < 2mb and Br(/ii — > aiOi) large enough 
to escape the LEP limits on the Z + 2b channel. For such points, the smallest allowed 
values of ma^, fnh^ and m/^g are about 120 GeV, 105 GeV and 180 GeV. We see that at 
low 772(12 < 200 GeV, ~ "'^102; while at large > 600 GeV one finds ~ fnh^. At 
intermediate rria^ , one can have either degeneracy. On the right-hand side of Fig. [211 we plot 
the relative coupling strength Cd of each of the heavier Higgs bosons as a function of its mass. 
Correlating with the mass plots, we observe that for rria^ < 200 GeV it is always the h2 that 
is degenerate with the 02 and that both have Cd ~ tan /3. For ^ 600 GeV, we always 
have m^g ~ and Cd{h-i) and Cd{a2) are both ~ tan/3. For in the intermediate mass 
range, the situation is more complicated and we can get cases where ~ rnh^ ~ m^g and 
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while Cd{a2) ~ tan/3 one finds that the /i2 and hs can share the tan/5 enhancement factor. 
This is illustrated in Fig. [251 

Overall, it is clear that if tan/3 is large a search at the LHC in channels such as gg 
bb + Higgs could reveal a signal so long as m^^ is not too large. For nia^ ^ 250 GeV, the 
Tevatron might also be able to detect this kind of signal when tan (3 is large enough. 



IV. CONCLUSIONS 



Eliminating EWSB fine tuning in supersymmetric models has become an important issue. 
In the NMSSM, this is most easily achieved by allowing the lightest Higgs boson, hi, to have 
mass of order 100 GeV, as naturally predicted after radiative corrections for stop masses in 
the range of a few hundred GeV. The modest stop masses imply no significant fine-tuning. 
However, such an hi is typically Standard Model like in its ZZ coupling and can escape 
LEP limits only if the dominant decay is hi — > aiOi (so that the hi bb decay is sufficiently 
suppressed to escape LEP limits) and if m^^ < 2mb (so that the aiai final state does not 
feed into the Z + b's final state that is strongly constrained by LEP data). In this paper, we 
have considered the degree to which the GUT-scale soft-SUSY-breaking parameters must 
be tuned in order to have m^^ < 2mfc and Br(/ii aiOi) > 0.7 (the rough requirement 
for suppressing the hi — > bb mode sufficiently). We have found that such a scenario need 
not have significant tuning. We began by assessing the tuning required of the m^-scale 
parameters Ax{mz) and A^^mz) that primarily control both and Br(/ii — * aiOi). This 
tuning was quantified via 

d log Tn^ 

FMAx = max{\FA,l\FAj} , = , (39) 

evaluated at scale mz- We found that so long as fi is not too large (|/i| < 200 GeV) 
then the values of A^{mz) and Axijnz) need only be tuned to a level of order Fmax ~ 
10 — 20, corresponding to tuning in the range of 5% to 10%, for the magnitudes of Af,{mz) 
and A\{mz) that are of order those automatically generated by radiative evolution from 
small GUT scale values. Further, these same RGE-generated values automatically give the 
required large values of Br(/ii aiai). We also discussed how these tuning estimates based 
on Fmax will generically greatly overestimate the tuning with respect to GUT-scale soft- 
SUSY-breaking parameters. In any SUSY scenario in which Ax{mz) and A^{mz) deriving 
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from RG evolution are dominated by a GUT-scale parameter or model-correlated set of 
GUT-scale parameters, generically denoted by p, the tuning with respect to p is given by 
fp ~ -^Aa + ^iid the Ax{Tnz) and Af^{mz) regions with < and Br(/ii — >■ 
Oiai) > 0.7 one finds F^^ ~ As a result /p can easily be very modest in size, even for 

quite large /i — the opposite-sign correlation becomes increasingly effective as // increases. 

Thus, we have demonstrated that GUT-scale parameters can be chosen so that fine-tuning 
can essentially be eliminated for the NMSSM scenario of a light Higgs with uih^ ~ 100 GeV, 
decaying primarily via hi aiOi with < and Br(/ii — > aiOi) > 0.7 (both being 
required to escape LEP limits). We regard this scenario (or something similar) as a highly 
attractive possibility for the Higgs sector. It solves both the /i problem and the fine-tuning 
problem. It does, however, introduce the need for proving viability of LHC discovery modes 
involving hi oiOi — > 4r or 4 jets, where the jets could be c, (7, s, . . .. The 4 jets mode 
is certain to be quite difficult to probe at the LHC - only diffractive pp — * ppX might 
provide a signal in the form of a bump in the Mx distribution at rrth^. In this regard, it is 
important to note that our work shows that minimizing the GUT-scale fine-tuning measure 
fp appears to mildly prefer m^^ > 2m^. Thus, a strong effort should be made to develop 
hi ^ flifli 4r discovery modes at the LHC and to see if LEP data can constrain this 
possibility. Of course, discovery of the hi at the ILC through a bump at Mx = nih^ in the 
e"'"e~ ^ Z* ^ ZX channel will be extremely easy. 
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